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Abstract 

We show that, under reasonable hypotheses, it is possible to study two-body non-leptonic 
weak decays in numerical simulations of lattice QCD. By assuming that final-state interactions 
are dominated by the nearby resonances and that the couplings of the resonances to the final 
particles are smooth functions of the external momenta, it is possible indeed to overcome 
the difficulties imposed by the Maiani- Testa no-go theorem and to extract the weak decay 
amplitudes, including their phases. Under the same assumptions, results can be obtained also 
for time-like form factors and quasi-elastic processes. 
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Introduction 



Exclusive non-leptonic decays are a fundamental source of information on quark weak inte- 
ractions and on the strong interaction dynamics. Unfortunately, a theoretical description of 
exclusive decays based on the fundamental theory is not possible yet. Over the years, several 
methods have been introduced to estimate the relevant matrix elements: vacuum saturation, 
bag models, quark models, QCD sum rules, l/Nc expansion, chiral Lagrangians, factorization, 
etc. None of these approaches is fully satisfactory. In the first three cases, the methods are 
in no way a systematic expansion of the fundamental theory in a small parameter. Thus the 
result of some calculation can be either correct or wrong, but it is impossible to reduce the 
uncertainty of the theoretical prediction by computing higher orders in the expansion. With 
QCD sum rules it is very difficult to improve the results in a systematic way because only few 
perturbative terms and matrix elements of higher dimensional operators can be estimated; 
in the expansion nobody has succeeded in computing consistently the matrix elements 
of the effective Hamiltonian beyond the lowest order and chiral Lagrangians have a rather 
limited domain of application. It may well be true that factorization holds for weak decays of 
hadrons containing a heavy quark, although a real proof that this is the case is still lacking 
and moreover there are amplitudes that cannot be factorized. 

The lattice approach has been used to obtain results based on first principles for a wide set 
of relevant physical quantities such as the hadron spectrum, the meson decay constants, the 
form factors entering in semileptonic and radiative decays, the kaon i?-parameter Bk, etc. For 
exclusive non-leptonic decays, however, this method is no better than the above-mentioned 
approaches, because of the Maiani-Testa No-Go Theorem (MTNGT) and the activity in this 
field H, |2| has completely stopped after the publication of ref . . 

In this paper we show that, under quite reasonable physical hypotheses, it is possible to 
extract predictions for the relevant matrix elements in numerical simulations of lattice QCD, 
in spite of the MTNGT. Our proposal has a wide domain of applications, ranging from kaon 
decays, K — > vrvr, to charm and beauty decays, D Kir, B Dtt, B vrvr, etc. If 
successful, it will also allow a check of several interesting theoretical ideas as the factorization 
of the amplitudes and the scaling laws that can be formulated in the Heavy Quark Effective 
Theory (HQET) when the mass of the quark becomes large. Under the same assumptions, 
it is possible to try a numerical calculation of the time-like form factors, e.g. those relevant 
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in e"'"e~ vr+vr^ or nn, and of those entering in quasi-elastic electron-nucleon scattering 
processes, e.g. e~ + p ^ e~ + A + vr. 

The MTNGT states essentially the following: 

• In the calculation of a two- (many-) body decay amplitude performed in the Euclidean 
space-time, which is the only possibility in Monte Carlo simulations, there is no distinc- 
tion between in- and out-states. As a consequence, the matrix elements that one is able 
to extract are real numbers resulting from the average of the two cases, e.g. 

{ttttIHwIK) = ^[^n{7T7r\nw\K) +out {ttttIHwIK)) . (1) 

Any trace of the phase due to final-state interactions is then lost. This jeopardizes the 
possibility of any realistic prediction for the matrix elements. For example, we know 
from the measured Ax/2 ^-iid A3/2 amplitudes in D — > Kt: decays that there is a phase 
difference of about 80°. 

• Matrix elements are extracted on the lattice by studying the time behaviour of appro- 
priate correlation functions at large time distances. Maiani and Testa showed that what 
can be really isolated in this limit are the off-shell form factors corresponding to the 
final particles at rest, e.g. {'k{ptt = 0)7r(p7r = ^)\'Hw\K). For kaon decays, we can use 
the chiral theory to extrapolate the form factor to the physical point[|. This is certainly 
not the case for D- and S-meson decays. In the latter case it is not possible to obtain a 
realistic prediction for the matrix element. 

We will show that both the difficulties raised by Maiani and Testa can be overcome under 
the hypothesis that final-state interactions are dominated by nearby resonances, the couplings 
of which to the final-state particles satisfy some smoothness condition]^. By varying the spatial 
momentum of the initial and final hadrons, it is then possible to extract the physical matrix 
elements, including the phase due to the strong-interaction rescattering of the final states. 
This allows also the calculation of the relevant parameters of the resonances. 

Two observations are necessary at this point. 
^ In the chiral limit final-state interactions are negligible and the problem discussed before is solved. This is thus 

the only situation in which the matrix element can be evaluated. 

^ A different approach, based on the study of the two-particle energy spectrum in a periodic box can be found 

in ref. g 
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The assumption that final-state interactions are dominated by resonances is not new and 
has been successfully applied to phenomenological studies of charmed meson decays, see for 
example ref. Q (together with the factorization hypothesis that we need not assume). 

At this stage, we do not know if the procedure that we are proposing can lead to useful 
results in practice. It remains to be seen whether, with a reasonable size of the lattice and 
number of gauge field configurations, it is possible to extract the matrix elements with a 
satisfactory accuracy. A feasibility study is currently under way on the APE machine. 



The three-point EucHdean correlation function 

Following ref. (|, we first examine the Euclidean three-point function G{ti,t2,q, —q) 

G{ti,t2,q, -q) = (o|r [ii^{ti)Ji_^{t2)H{Q)\ |o) = (o|n,-(ti)n_,-.(t2)i^(o)|o) , (2) 

when ti> t2> 0. In eq. (^), ^q{t) is an interpolating field of the final-state particle (denoted 
as "pion" in the following) with a fixed spatial momentum 

Ii^{t) = [ d^x e-^«-%(f , t) ; (3) 



i^(0) = H{x = 0, t = 0) is any local operator that couples to the two pions in the final state; 
T [. . .] represents the T-product of the fields and the vacuum expectation value corresponds, 
in a numerical simulation, to the average over the gauge field configurations. 
When ti ^ CO 

G{tiM,q,-q)^Y.^m?{ti)\n){n\Ii_^{t2)HmQ) ~ ^e-^?*^G3(t2) , (4) 

where 



= + , = (0|n(0)|g-) (5) 

and 

G^{t) = {q\ii-<im{m (6) 

for t > 0. 

Inserting a complete set of out-states, we can write 

{q\U-g<t)HmO) = Y,{2Trf6HPn) {q \e'''n{0)e-^'\n) ^ut out{n\H {0)\0) = 
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5:(27r)3<53(p„)(g-|n(0)|n)o„io.i(n|i/(0)|0)e-(^"-^^)* 



(7) 



2E^ 

J2{27TfSHPn){mmn)outout{n\HmO)e-^^--^'r'^' 

- ^ J connected 

where the term proportional to out{Q, — ^-f^(0)|0) on the r.h.s. is the disconnected contribution. 
We now make the following assumptions: 

1. The only possible eigenstates of the Hamiltonian are n-pion states. 

2. The interaction between pions is dominated by a narrow resonance (denoted as a in the 
following) exchanged in the s-channel. 

3. The coupling of a to the pions is a smooth function of the external momenta. 

Under hypotheses 1. and 2., the only possible intermediate states are two-pion states. We 
thus obtain 



Gs{t) = Y,i'^TTf6^{Pn){q\Umn)outout{n\H{0)\0)e 



-{En —Eg)t 



2E. 



+ 



^(27r)3<53(P„)(g|n(0)|n)o„t out{n\H{0)\0)e' 



{En — E^)t 



(8) 



connected 



In eq. (P) = 2E^ ; \n)out = \k,—k)out] y/s = E = 2£'g; the sum over the intermediate 
states is given by|^ 



J2{27Tf6HPn 



dE 



2tt J (27r)32S^^ (2^)32^^^ 



. ^, . (2TrY6-\ki + k2) = 

(2^)32i?,^^ (27r)32i?,^/ ^ ^ ' 

{27r)^6'\ki + k2)S{E^^+E^^-E) 



d^ki d^k2 /^_x4.3 



and we have used 



out{q,-g]HmO) 

We now define 

{q\U{0)\k,-k)out = 



M^-Sg-iX{Sgy'''=P«- 

2VZ^ \V{s)\^ 



(9) 



(10) 



{E + 2E^-ie){-E + 2E^-ie) - s - iX{s) 

\V{s)\-' 

{E + 2E^ + ie){-E + 2E^ + ie) - s + iX{s) 



(11) 



further factor 1/2 in the case of identical particles is understood. 



In eq. ( pi] ) we have introduced the a-vr-Tr couphng V{s). As stated before, we assume here 
that the couphng is such a smooth function of the external momenta that we can use the 
same "physical" coupling also for the off-shell pion, which is annihilated by 11(0) in the matrix 
element |n(0)|/c, —k)out- In general we should write 
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out{k,-k\U{0)\q) 
\V{s)\' 



M2 



iX(s) 



= VZu X 
X ^( 1 - 



{E + 2E^-ie){-E + 2E^ 
E \ 



2E. 



(12) 



with the condition that the modulating factor of the off-shellness 1 — E/2E^ satisfies the 
condition ^(0) = 1 (see below). The factor in parenthesis is (up to a factor 1/2E^) the 
propagator of two non- interacting pions [0 . 



With the definition in eq. (11), the quantity in square brackets satisfies 

2^ \V{s)\'' 



lim {-p^ + Ml) 



M2 



{E + 2E^-ie){-E + 2E^-ie) - s - iX{s] 



, (13) 



which is the usual LSZ-reduction formula for on-shell particles, cf. eq. (15) of ref. |^. In 
eq. (|13|), p2 = E'^ + — 2E Eg' is the squared four-momentum of the off-shell pion. Combining 
eqs. (^ and ( pi]) and using the identity 

1 



E - 2E^ - ie 



V 



1 

E - 2Efi 



+ mSiE - 2E^) 



qj 1 



(14) 



we obtain 



Zn 

-e 9 X 



2E, 



aisq) 



M2 



dE 



^e^9^* X 



Sq - iX{Sqj 

9{s) 



Xis) 



-Et 



2E, 



^1 



{E + 2E^ + ie){-E + 2E^ + ie) [W^ - sf + X'^ {s) 

( ajSq) ig{Sq)X{Sq) \ ^^E^t 

I M2 -Sq- iX{Sq) {M^ - Sq)^ + X^{Sq) i 



\ f+oo 



dEe 



-Et 



9{s)p{s) 
AEl - E^ 



(15) 



where, in the integrand of the r.h.s., we have denoted as E'^ or s the same quantity, s = E"^. 
This will facilitate the comparison of the above formulae with the more general case of a 
two-pion final state with total momentum different from zero, cf. eqs. (Eol) and (ETI) below. In 



eq. (15) we have used the relation X{s) = Mo-r(s), where r(s) is the "width" of the cr-meson 



defined as 

The spectral function p{E'^) is given by 



{M^-E^y + x{E^y ' 

which becomes 

1 T{E^)/2 

2M^ (M, - + (r(£;2)/2) 

for a narrow resonance. We also have 



piE') = — . 1.^.. ,..2 (17) 



2M, 
as r ^ 0. 

We then establish the relation 

a{sq) ig{sq)X{s 



p{E^)^^S{M^-E) (18) 



M2 -Sg- iX{Sg) (M2 - Sg)^ + ^^(s,) 

9{Sq) ^ 1 ^(Sg) giSg) 



M^-Sg- iX{Sg) 2\M^-Sg + iX{Sg) - S g - iX{Sg) J 

out{q,-mmO) + l(^n{q,-q\HiOm " out{q,-mmO)) , (19) 



2 

in agreement with eq. (21) of ref. Q. The last term in eq. (|l5|) corresponds to the last term of 
eq. (23) of ref. |^]. The definitions of out{Q, —q\H{0)\0) and in{Q, — ^-f^(0)|0) are in agreement 
with the standard relation between the strong rescattering phase and the parameters of the 
resonance a 

where the phase 6{s) is defined up to an ambiguity of vr = 180°. We choose < 5{s) < it. 

In the case of non-interacting pions, out{Q, ~Q\H{0)\0) = in{Q, —q\H{0)\0) = y/Z^, and one 
finds 



G{h,t2,q,-q) = i^j e-^^^(*^+*^)yz;. (21) 



Equation (|2l|) has a very simple interpretation. The two non-interacting particles are created 
in the origin t = with amplitude ^/Z^ and propagate from the origin to t2] from t2 to ti 
only one particle state propagates. 

In the limit of a zero- width resonance {X{s) 0) we obtain 



Giti,t2,q,-q) 



2E, 



(22) 



which is in agreement with the result, obtained for a narrow resonance, given in eq. (3.3) of 
ref. Q. We notice that eq. (3.3) of ref. Q is really valid only at small time-distances. According 
to the MTNGT instead, for any non-zero value of the width, the correlation function at large 
time-distances is always dominated by the state with the two pions at rest. 

For further use, it is convenient to introduce the quantity R^{t2,q} defined from the 
relation 

„ 

(23) 



Gih,t2,q,-q)=(^^^ e-^^<'^^'^^R^{t2,q) 



where the label H in R^{t2,q) denotes the local operator used to create the two pions. 

It is straightforward at this point to derive the expression that enters in a two-body non- 
leptonic decay of a meson (which we will denote as D in the following). The starting point is 
the four-point correlation function 

G{ti,t2,tD,q,-q,PD = 0) = ^ lim {Ilg{ti)U_g-{t2)nw{0)Di Md)) 

= (^)'e-^^^(*^+*^)i?^-(t2,g1^e-^^-l*-l, (24) 



with 



R^^{t2,q) 



+ 



1 



M^-Sg-iX{s,) 2\M^-Sg + iX{s,) M^-Sg-iXisg) 



+ 



and 



vr AEl -E^ 2i\M^-s- iX{s) - s + iX{s) 

out{k,-k\nwmD{pD = 0)) 



(25) 



(26) 



M2 - s - iX{s) ' 

where Tiw is the weak Hamiltonian. We have not considered complications coming from 
the presence of (axial-)vector particles and different flavours. The formulae reported in this 
section can easily be extended to these cases. 



7 



The Watson theorem and physical interpretation of 
the results 



Let us start from 



G{h,t2,q,-q) 



9[Sq) 



2E. 



9{Sq) 



9{Sq 



2\M^-Sg+ iX{Sg) M^-Sg- iX{Sg) 



2M^ TT 



M^-Sg- iX{Sg) 
+ 

9{s) 



(27) 



9{s) 



4Ej-E^2i\M^-s-iX{s) M^-s + iX{ 



and introduce the following relation 

out 

{k,-k\Hmo) 



9is) 



(28) 



M^-s- iX{s) 

The Watson theorem ensures us that (in the absence of CP violation) A(s) is a real quantity 
0. It follows that 



g{s) = JiM^ - s + iX(s))(M2 - s - iXis)) x A{s) 



is also a real quantity. Equation (pO[) implies 



We may then write 



.(^V^>(0)|0) = ^,_f;\^^^^ =e-^(-)A(.). 



(29) 



(30) 



(31) 



with 



A{sg)cos6isg) + {AE^)V 



cos (5(s) 
sin (5(s) 



2M„ ^ 



AEl - E^ 



A{s) sin 6{s) 



Ml 



^{M'^-sY+X{sf 
X(s) 



^{M'^-sY+X{sY 
Notice that at threshold 6{s)\g^^]yp = 0, cf. eq. ([l^). 



(32) 
(33) 
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Following ref. we can compute the leading behaviour of G{ti,t2,q = 0,—q = 0) for 
large t2 and express the result in terms of the parameters of the resonance. We obtain 



G{ti,t2,q = 0,-q = 0) 



ami' 



+*2)^(4M2) 



\ — a\ 



Trt2 



where the scattering length a is given by 

^ _ |y(4M2)|2 



levrM^ M2 - 4M2 ' 



(34) 



(35) 



in agreement with eq. (10) of ref. Q^. 



Applications to lattice calculations 



For Z)-decays eq. (31) becomes 

G{ti,t2,tD,q,-q,PD = 0) = 
\A^{sg)cos6isg) + {AE^)V 



2E. 



e 1 



dE_^_f^E_2E^)t2 1 
2M^ TT 



2Md 



MD\tD\ 



(36) 



where {s) = V(M^ -s + iX{s)){M'^ -s- iX{s)) x A^{s). 

For numerical applications, it is convenient to consider the amputated correlation function 
given by the ratio 



R^'^{t2,q) 



G{ti,t2,tD,q, -q,PD = 0) 
Sn{ti, E^) Su {t2 , E^) Sd {to , Md ) 



where 



Sj,{t,,E^) = ^e-^^'^, 



2E. 



(37) 
(38) 



and similarly for the other meson propagators; 
R^^{t2,q) = \ (sg) cos 6{sg) + {4E^)r 



°° dE ,E-2E^)t2 1 



2M^ vr 



AEl - E^ 



A^{s) sin 5{s) 



(39) 



Until now, for simplicity, we had chosen = 0; the generalization of the above formulae 
to pd 7^ 0, corresponding to G{ti,t2,tD,qi, —q2,PD = qi — q2), is straightforward: 



R^'^{t2,qu-q2) 



G{ti,t2,tD,qi, -q2,PD = qi- q2) 
•Sn (i 1 , Eff^ ) Sji {t2 , Eff^ )SD{tD, E.^^ ) 



There is a difference of a factor of two, due to the fact that in our case the particles are distinguishable. 
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(st) COS 6 (st) + {2Et)V 



dE_ ,E-ET)t2 1 aW 



E TT 



A'^{s)s\n5{s) 



where 



Et = E^^ + % ; Emin = + (qi - Q2)2 ; st = E't- (qi - q2) 

and s = E'^ - {qi- q2)'^. 



(40) 



(41) 



We now explain the strategy to extract the physical information, i.e. the matrix element 
of Hwi including the phase, from R^^^' {t2,q)- The extension to i?^^(t2, ^i, — 92) is straight- 
forward. 

Let us imagine to make a calculation on a lattice of volume x T. The lattice version of 
eq. (|l5|) is given by|: 



R^^{t2,q)= A{sg)cos5{sq) + 



4E„ 



TT 



E 



1 



AEl - Ef 

q I 



A{s) sin 6{s) 



(42) 



where all the quantities are given in units of the lattice spacing, 



Ei 



27r 



^s = 2E^ with k = —{n^,ny,nz); 



(43) 



Sq has been defined before; nx,y^z = 0, 1,... ,L — 1 and Y^'e^ denotes the sum over all the 
values of the energy corresponding to the momenta k allowed by the discretization of the 
space-time on a finite volume, excluding those corresponding to Ei = 2E^. Different combi- 
nations of momenta corresponding to the same energy should be included only once in the 



sum appearing in eq. (42), since the factors multiplying already account for the phase 
space integration, at fixed total energy. Thus, for example, \k\ = 2tt/L corresponds to six 
possibilities^, (±1,0,0), (0, ±1,0), and (0,0, ±1), but has to be counted as a single term in 
the sum. AE'j = EiJ^i — Ei is the difference of the nearest successive allowed values of Ei 
{Eq = 2M^, El = 2y/M'^ + {2iT/Lay^, E2 = 2^/M'^ + 2(27r/La)^ etc.). One can show that 
the expression in eq. ( p2[ ) tends to the corresponding continuum one in eq. ( |3^ ) as L ^ 00. 



^ We assume here that the discretized version of the two-body phase space is a good approximation of the 
continuum one. In this case we can change the integrals into sums over discrete values of momenta or energies 
without problems. Further technical complications arising from the different density of states in the continuum and 
on the lattices currently used in numerical simulations will be discussed elsewhere. 

^ In the following, k = 2tt/ L{nx,ny,nz) will be simply denoted by (1,0,0). 
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By varying the external momenta we can study 



tanS{sg) = ^ and \A{sg)\ = ^ Al + (44) 



as a function of the centre-of-mass energy, where Ac = A{sq) cos 5{sq) and As = A{sq) sin 5{sq). 
From the behaviour of tan5(sg) as a function of Sg, we can reconstruct the mass and width 
of the resonance and extrapolate the amplitude to the physical point. 

Ac and As can be found by studying the dependence of R^^' {t2,q) on t2, at fixed q, for 
ti > t2 > 0. 

Several observations are important here: 

• The range of values of Sq, which is available in current numerical studies, is limited 
because the momenta on the lattice are quantized and, in order to avoid discretization 
errors, \q\ ^ 1 {\q\a <^ 1, where a is the lattice spacing, if we do not work in lattice 
units). In order to enlarge the range of values of Sq, it is convenient to have a D-meson 
with non-zero momentum, corresponding to qi ^ q2- 

• To derive the properties of the resonance, it is not necessary to use the four-point correla- 
tion function. We can instead study the three-point i^-II-n correlation G{ti,t2, qi, — ^2) 
as a function of the momenta. There are several advantages in doing so. First of all, 
an additional information on the parameters of the resonance can be exploited when 
analysing the four-point correlation function. Moreover the signal of the H-H-U corre- 
lation is expected to be much less noisy than in the four-point case. A final advantage is 
that, by varying the quantum numbers of the field H{0) we can excite and study different 
resonances, for example a vector-like one. 

We now present an example to show how [t2, q) would appear as a function of ^2 for 
specific values of the mass and width of the resonance. What happens is the following. In the 
absence of final-state interactions, R^^{t2,q) would display a plateau as a function of ^2- In 



the presence of the resonance, the second term in eq. (42) increases (decreases) exponentially 



in t2 for states with Ei < 2E^ [Ei > 2E^). It can easily be seen that only the states with 



Ei < 2E^ give an appreciable contribution: the last term on the r.h.s. of eq. (42) gives 
a dramatic effect on the plateau for large values of t2- This is shown in fig. 1, where the 
following parameters have been chosen: 

1. L = 24 and T = 64; 
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2. f3 = 6.2, corresponding to a ~ 2.9 GeV; 

3. a resonance with a mass = 1.9 GeV and a width r(M^) = 0.3 GeV, as was used for 
D-decays in ref. 

4. r(.) = r{M^)^[{s - 4M^)M^]/[(M^ - Ur^)s]; 

5. the masses of the pseudoscalar mesons have been taken to be M-,^ = 0.4 GeV; 

6. the external momenta are qi = (1, 1,0) and q2 = (—1, —1,0), so that there are only two 
intermediate states, of zero total momentum, corresponding to A; = and k = 2tt/L, and 
energy less than the energy of the external state; of the two states, the one with the two 
pions at rest does not contribute since 5(4M^) = 0. 

7. (s) has been taken to be constant and equal to 1. 

By fitting the curve of fig. 1 as a function of ^2 it is possible to extract A{sq) cos{sq), with 
Sq = 4[M2 + 2(27r/La)2], and ^(sfe) sin(sfc), with Sk = ^[M^ + {2'K/Laf]. In the same way, 
by increasing the energy of the external two-pion state, and using also states with non-zero 
momentum, we can extract Ac and Ag as a function of Sq and extrapolate the amplitude and 
the phase to their physical value. 

Conclusion 

We have shown that, in spite of the MTNGT, it is possible to extract the relevant information 
on two-body decay matrix elements, under the hypotheses that final-state interactions are 
dominated by nearby resonances and that the couplings are smooth in the momenta. We have 
also outlined the strategy to extract this information from the Euclidean correlation functions 
that can be computed in standard numerical simulations. The feasibility of our proposal 
requires a dedicated study, which will be presented elsewhere. 
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Figure 1: X{t2) = R^'^{t2,q) as a function of t2 with the parameters given in the text. AEa — 
2{E^- E^)a, where 2E^^ ^ = 2^M2 + 2{2t: /Laf and 2E^^^^ 2^ Ml + {^-k j Laf . 
and are the corresponding phases. The dot-dashed line is the constant term corresponding to 
{sg) cos6{sg) , the dashed curve is the contribution of (sk) sin d{sk) and the solid curve is the 
sum of the two terms. 
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